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Abstract 



We comment on recently discovered 1/Nc corrections to the nucleon axial current in the 
framework of the Nambu-Jona-Lasinio soliton. This kind of corrections arises only for 
a special treatment of ordering ambiguities of collective coordinates and operators in the 
semiclassical quantization. In addition to the missing derivation of this special quantization 
scheme from first principles its naive application violates partial conservation of the axial 
current (PCAC). We show how within this scheme PCAC can be restored and determine 
the corresponding I/Nq corrections to the equation of motion for the chiral soliton. The 
resulting self-consistent solution allows to evaluate the nucleon axial coupling constant qa 
directly as a matrix element of the axial current as well as indirectly from the pion profile 
function. Enhancement is found for those baryon properties which are sensitive to the long 
range behavior of the pion profile. 
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Chiral symmetry and its spontaneous breaking are generally accepted as guiding principles 
for the construction of models describing hadron physics. A prominent feature of chiral 
symmetry is the partial conservation of the axial vector current (PCAC). PCAC relates the 
divergence of the axial vector current A"^ to the pseudoscalar-isovector pion field tt": 

d'^A; = mlUn\ (1) 

Here m.,^ = 138MeV and f-^ = 93MeV denote the pion mass and decay constant, respectively. 
Since the axial current is of the nature of a Noether current its (partial) conservation heavily 
relies on the pion fields satisfying their equation of motion; i.e. eqn. (|I|) is equivalent to 
the equation of motion for the pion field. Therefore variations of the axial current imply 
corresponding changes in the equation of motion for the pion field. 

In baryon physics the axial current plays an important role since its matrix element 
between nucleon states, the nucleon axial coupling constant qa, directly relates to the neutron 
(3 decay amplitude. It has been a long standing problem for chiral soliton models of baryons 
that the predicted value for qa is only of the order 0.6 to 0.9 underestimating the experimental 
value 1.26|l|0 Recently it has been claimed that corrections might provide the 

missing amount to match the empirical value. However, these corrections are only present if 
a certain ordering prescription for the generators appearing in the semiclassical quantization 
scheme is adopted. We will resolve this point below. Additionally, the way these corrections 
are treated in ref. |Q violates PCAC. In that calculation only the corrections of the axial 
current were taken into account disregarding the corresponding changes of the soliton profile. 
In this letter we will demonstrate that these corrections are mitigated but nevertheless sizable 
when PCAC is properly accounted for. 

Let us briefly review the connection between the equation of motion for the meson profile 
function and PCAC in chiral soliton models. The static solution to the Euler-Lagrange 
equations acquires the hedgehog shape for the non-linear representation of the meson fields: 

?7o(r)=exp(zT-f e(r)), (2) 

0(r) being the chiral angle. The equation of motion for 0(r) may generally be expressed in 
the form ^(0) = 0. The hedgehog configuration (|^) is quantized by introducing collective 
coordinates for the (iso-) rotational zero modes: 

U{r,t) = R{t)Uo{r)R\t), R{t) G SU{2). (3) 

Due to the vanishing time component of the axial current for the configuration (^) PCAC 
reduces toB 

d,A'^ = DathH^) - f^mlDabhsinQ. (4) 

Dab = (l/2)tr(Tai?Tbi?"'') deuotes the adjoint representation of the collective rotation. The 
PCAC relation is obviously recovered if the chiral angle satisfies the appropriate equation of 
motion which is obtained upon inversion of eqn. (^): 

= ^(0) = hDab^^A'^ + />'sine. (5) 
"Note, however, that in the chiral quark model gA is somewhat overestimated [||. 

^For our sign conventions the pion field of the hedgehog configuration is given by tt" = — Dab ftsinO- 
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The axial coupling constant gA is defined as the matrix element of the axial current 
between proton states at zero momentum transfer: 



gA = 2{pUAl\p^). (6) 

The defining equation for the pion decay constant provides the asymptotic form of the axial 
current: 



hm At = U lim <9,7r'^ = -fl lim (D^khO) ■ (7) 



The identity dj{riA'^) = A" + TjOiA^ then allows to determine qa from the pion profile via 
PCAC m): 



47r „,ae 



3''^iJ-oo dr 



Stt 



9A = —f: lim R'^ _^ - —f:mi / rfrr^sinG (8) 



--R 9 



wherein extensive use has been made of the matrix element (-D33)proton = ~"l/3 in the two 
flavor case. Note the appearance of the factor 3/2 in the surface term which arises at zero 
momentum transfer in the chiral limit (m7r = 0) We would like to emphasize that the 
equivalence of eqns. and (|^) manifests a direct consequence of PCAC for chiral soliton 
models. Phrased otherwise: PCAC allows to read off qa directly and uniquely from the pion 
profile. 

So far these considerations have been fairly general and apply to various kinds of chiral 
soliton models. Let us now turn to the specific case of a two flavor Nambu-Jona-Lasinio 
(NJL) model of scalar and pseudoscalar interactions in the isospin limit [^]. After bosoniza- 
tion the action Anjl may be expressed as a sum Anjl = Af + Am of a fermion 
determinanto 

Af = Tr log^ (ip) = Tr log^^ {i^ - Mt/^^ ) (9) 
and a purely mesonic part 

^™ = ^ / d'xti (f/ + t/t - 2) (10) 



where G is the NJL coupling constant. For a finite pion mass it may be eliminated via the 
relation G = mM/m'^f'^. For a given constituent quark mass M the cut-off A is determined 
such as to reproduce the empirical value of the pion decay constant while the current quark 
mass m is obtained by solving the gap equation 0. Throughout this letter we will apply 
Schwinger's proper time regularization prescription to the diverging fermion determinant. 
The NJL model supports static solitons (0) which are obtained as the solutions of the Euler- 
Lagrange equations to the static energy functional p, ^ |T0|. 

Only the fermion determinant (||) contributes to the axial current. Since Af = ^vai + ^^sea 
splits into valence and sea quark parts we have a corresponding decomposition of the axial 
current: 

j^a ^ ^(val)a ^ ^(sea)a_ ^^^^ 



'^For simplicity we treat only the chiral field U as space-dependent quantity and keep the scalar field equal 
to its vacuum expectation value (</>) = M . 
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The valence quark part is given by 

= iVc^(.al)7.75 Y^(val) (12) 

with Nc being the number of color degrees of freedom. In order to determine the valence 
quark wave function \l/(vai) we need to construct the eigenfunctions ip^^ of the static Hamil- 
tonian /iq: 

hoiP^ = {oi-p + fSMU-^') = e^iP^ . (13) 

Denoting the corresponding eigenstate and the state with the lowest positive energy 
eigenvalue |0), the valence quark wave function \l/(vai) is given as a perturbation expansion 
in the angular velocities {i/'2)T • O = K''{t)R{t): 

v^M) = rL + Ij: + OiO^). (14) 

Substitution of eqn. (|14[) into eqn. (0) yields 

^ ^ Mt^o ^0 - ^ 

+ hD^,n,)J:^^^4<y^l4'^, + 0{f2'). (15) 

The terms linear in f2 on the RHS of eqn. (^) give vanishing contribution to qa as long as 
the angular velocities are considered to be pure c-numbers; i. e. there is no I/Nq correction 
to gA if is treated as a classical quantity. However, this is no longer the case when the 
semiclassical quantization prescription 

is imposed. Here Ji are the spin operators in the space of the collective coordinates and 
denotes the (iso-) rotational moment of inertia. The analytic expression for in the NJL 



model may be found in ref. |TT|. Using furthermore the commutation relation 

[Jj, Dab] = i^jbmDam (17) 

the collective operator involved in the correction to qa is then found to be of the same 
structure as in the leading term in agreement with more general considerations [|r^. Up to 



linear order in the valence quark contribution to the nucleon axial coupling constant is 
obtained to be 0: 



As is of the order A'',^ the correction in (|T8D is of order 1/Nc- The result ([T8|) is, of course, 
ambiguous. This is due to ordering ambiguities of the semiclassical quantization prescrip- 
tion([T6|) arising when the classical collective coordinates are substituted by the collective 
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generators. Since there is no derivation of this prescription from first principles uncertain- 
ties of this kind are unavoidable. Ignoring, for the time being, this kind of ambiguities we 
will explain the relevance of PCAC once this special quantization prescription is adopted. 
In order to derive the corresponding 1/Nc corrections to the equation of motion we need to 
evaluate the divergence of (0). Using the eigenequation (|I3]) we find: 



N, 



-fbDabdiA. 



(val)a 



C 



Mi/jqIS (— sin6 + ir ■ f75Cos9) ipo 
+ ^E^^^l/?75[r.f,r,#oCose. 



(19) 



The first term on the RHS of eqn. (19) is just the valence quarks' contribution to the 
classical equation of motion ^ for chiral angle while the second term represents the I/Nq 
correction we are looking for. According to ref. no I/Nq corrections to qa arise from the 
sea-part of the action ^sea^ 



sea 

9a 



^''■Eerfc(|^|)(/x|S3r3|/i) 
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(20) 



and thus there is also no further change to equation of motion. Note that in the chiral limit 
(m^ = 0) an overall factor 3/2 appears in eqns. (|TB|,^Dp |Q. The absence of I/Nq corrections 
to g^^^ has been doubted and small corrections may indeed be present [0. We will not go 
into this point further but rather concentrate on the 1/Nc corrections to the equation of 
motion stemming from the valence quark wave function since this suffices to illuminate the 
relevance of PCAC. The complete equation of motion we would like to solve then reads 



cos9(r) 



\Tj\fx) 



sine(r) dn(tPoi^o + tT{pT)) - 



Nc M 



(21) 



wherein p^*^ denotes the sea contribution to the scalar density defined in ref. [^. 

As the moment of inertia is a functional of the meson profile the appearance of 
in (l2l| ) transforms the original algebraic equation into an integral equation. Its numerical 
solution is obtained by iteration. For a given the self-consistent solution to eqn s. (P^ and 
(^Tj) is constructed. This solution then yields a new to be substituted in eqn. (|21|) . Since 
the resulting deviates only moderately from the one obtained when the 1/Nc correction 
is omitted a converging solution is already obtained after a few iterations. 

As can be seen from fig. (1) the pion profile at large radii is increased if the 1/A^c correction 
is taken into account. This increases the moment of inertia which in turn decreases 
the correction term in the equation of motion (|2l|) and in the expression for qa (0)- This 
behavior guarantees the stability of the 1/Nc corrected soliton and limits the change of most 
quantities, see table (1). E.g. the soliton energy Etotai is increased by a few percent 
only. As the slope of the soliton profile close to the origin becomes moderately larger, the 
valence quark energy eo decreases accordingly. Also the moment of inertia increases by a 
small amount yielding an even smaller nucleon-A mass splitting Ma - M„ = 3/2a'^. A 
similar behavior is found for the isoscalar root mean square radius (r^)){^^Q. In the chiral 
limit it is even overestimated when the 1/Nc correction is taken into account. In contrast to 
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Fig. 1. The chiral angle as the self-consistent solution of eqns. (13 2^): solid lines. The dashed lines refer 
to the omission of the \jo? term in (21). Displayed are the cases m.,^ = (a) and = 138MeV (b). 



Table 1. Static properties for the soliton solution of eqns. (^ 21): '0 corr'. '0 not corr' refers to the 
omission of the \j 0? term in the equation of motion ( ]2l| ) . The superscript for gji^ denotes the order of 
Xjo? in eqns. ( [l8| ) and (po[). (?^°*'* = + is the final result for the nuclcon axial coupling con- 
stant, eo, Efotah and (r'^)j{^Q are the valence quark eigenenergy, total energy of the soliton, moment of 
inertia and the isoscalar root mean square radius, respectively. We have used the constituent quark mass 
M = 400MeV. 





m^r = 


= 138MeV 




not corr. 


corr. 


not corr. 


corr. 




0.83 


0.82 


0.75 


0.81 


si" 


0.47 


0.35 


0.39 


0.30 


„(tot) 

9a 


1.30 


1.17 


1.14 


1.12 


gA from eqn. {i) 


0.82 


1.13 


0.76 


1.14 


eo/MeV 


190 


151 


215 


167 


Etotal/MeY 


1203 


1233 


1253 


1266 


aV(l/GeV) 


6.49 


7.50 


5.86 


6.09 




0.90 


0.98 


0.75 


0.80 
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the moderate change of other quantities the first order correction for drops significantly. 
Whereas in the chiral hmit this change is reflected also in the total sum for there is 
almost no net effect in the case of the physical pion mass because the zeroth order value for 
gA increases. 

Table (1) reveals also that without correction for the meson profile the value of gA differs 
significantly from the value as obtained using eq. i.e. PCAC is strongly violated. As 
expected the PCAC relation @ provides values which are equal to the zeroth order con- 
tribution. On the other hand, with the corrected meson profile PCAC is fulfilled within 
numerical uncertainties. Note that the extraction of gA from (^) becomes the more difficult 
the larger the soliton is, as all numerical calculations are done in finite boxes. We would also 
like to mention that we find some discrepancies from the numerical results of ref. 0. Some 
of these discrepancies can be resolved by multiplying their result for g^2^ for m„ = by 3/2. 
The necessity of this factor has been pointed out in ref. 

In conclusion, we have calculated the NJL soliton including 1/Nc corrections as they 
result from the ambiguous quantization "recipe" (0). As claimed already in ref. these 
corrections solve the problem of a too small value of the axial coupling of the nucleon gA in 
hedgehog soliton models. Taking into account these corrections in a way that PCAC is still 
fulfilled decreases the value of the correction for gA- Nevertheless, the total result is less than 
ten percent off the experimental value. Here we argued on the basis of PCAC that the soliton 



profile has to change also if we adopt the quantization ([T6|JT7|) . It would be interesting to 
construct the energy functional which corresponds to that operator ordering and therefore 
provides the meson profile calculated in this letter from a variational principle. Since the 
soliton mass is of order Nc this correction has to be of order A^^. On the other hand, as long 
as the existence of such a formalism cannot be shown, the application of "recipe" (|16|,|1^ 
to different orderings of collective coordinates and operators will at least be doubtful. After 
all, it would be somewhat astonishing if the problem of a too small gA for hedgehog soliton 
models turned out to be only a result of a peculiar quantization. 
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Table caption 



T3|,|2T1): 'O corr'. '6 not corr' 



Table 1: Static properties for the soliton solution of eqns 
refers to the omission of the term in the equation of motion (0). The superscript for 
denotes the order of l/a^ in eqns. ([T8|) and 
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g^^^ + g^^^ is the final result for the 



nucleon axial coupling constant, eo, Etotai, o? and (r^)}{^^Q are the valence quark eigenenergy, 
total energy of the soliton, moment of inertia and the isoscalar root mean square radius, 
respectively. We have used the constituent quark mass M = 400MeV. 
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Figure caption 



Figure 1: The chiral angle as the self-consistent solution of eqns. ([T3|j2T|) : solid lines. The 
dashed lines refer to the omission of the l/a^ term in (|2T|). Displayed are the cases m^r = 
(a) and = 138MeV (b). 
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